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The second-order coherence of photons scattered from a 
trapped Bose-Einstein condensate is found to be enhanced for 
the scattering angles that are either the same or symmetri- 
cal with respect to the direction of laser propagation. The 
enhancement occurs periodically with respect to the time in- 
terval between photons, and becomes less and less pronounced 
with increasing the strength of the interaction between atoms. 
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I. INTRODUCTION 

Realization of Bose-Einstein condensation (BEC) in 
trapped atomic gases 0-^ has created an interdisci- 
plinary sub field of physics called coherent matter wave. 
Developments of optical imaging techniques such as ab- 
sorption imaging and a nondestructive in situ method ||] 
have enabled us to observe interference between matter 
waves Q and dynamics of BEC 0,^- In a remarkable 
recent experiment, Hau et al. demonstrated that it 
is possible to make photons propagate through BEC at 
a speed of 17 m/s, which suggests that the index of re- 
fraction of BEC is 14 orders of magnitude greater than 
that of a glass fiber. In the present paper we predict 
yet another unique feature of BEC concerning correlated 
photon scattering. 

Various aspects of the optical response of BEC have 
been predicted p0|-|20|: a broad resonance |13 1^, an ex- 
tra structure in the spectrum |p^ , p^ , response to intense 
short pulses ||l^,|l^, non-Lorentzian line shapes p^ , pO[ , 
and a refractive index of a homogeneous gas ||l8| . In the 
above work, kinetic motion of atoms is either ignored or 
not explicitly taken into account because its time scale is 
usually much slower than that of optical processes. How- 
ever, optical signatures of BEC are expected to appear 
also in the time scale of kinetic motion, since the char- 
acteristic feature of BEC is the degeneracy of the kinetic 
degrees of freedom as well as the atomic internal degree 
of freedom. 

The aim of the present paper is to study correlations of 
photons scattered from trapped BEC that is irradiated 
by a weak far-off-resonant laser field, and to show that 
they exhibit some characteristic features associated with 
kinetic motion of the atoms in the trap. By scattering a 
photon, the condensate suffers a recoil, which affects the 
scattering of the subsequent photons. As a consequence 
a correlation arises between a pair of photons with re- 
spect to their scattering angles and time intervals. We 



show that the presence of the condensate enhances the 
probability of a pair of photons to be scattered into ei- 
ther the same direction or symmetrical directions about 
laser propagation, and the enhancement occurs periodi- 
cally with respect to the time interval, where the period 
is set by the frequency of the trap. 

The interaction between atoms cannot be ignored, for 
the profile of the condensate is significantly altered even 
for a few hundred atoms 22 1 . While in Ref . effects 



of the atom-atom interaction on the photon scattering are 
taken into account only through the atomic density pro- 
file, the dynamical effect of the interaction on the optical 
response has not been discussed. We take into account 
the atom-atom interaction using the Bogoliubov approx- 
imation ||2^ , and show that the correlations between the 
scattered photons are deteriorated by the interaction. 

This paper is organized as follows. Section |l| for- 
mulates the problem and derives an expression of the 
electric-field operator far from BEC irradiated by a laser 
field. Section III studies coherence properties of scattered 
photons. Section IV discusses effects of the atom-atom 
interaction on the scattering of photons, and Sec. ^ con- 
cludes this paper. 



II. FORMULATION OF THE PROBLEM 

We consider a situation in which trapped atoms are 
irradiated by a laser field which is sufficiently detuned 
from the resonant frequency of the atoms so that each 
photon is not scattered by the atoms more than once. 
Photons are scattered from the atoms and detected by 
photon counters. Quantum-statistical properties of the 
scattered photons can be studied through the electric- 
field operator E(r,t) in the Heisenberg representation. 
We derive here an expression of E(r, t) at the position far 
away from the trapped atoms compared with the spread 
of their wave functions. 

The Hamiltonian for the photon field is given by 



(1) 



where a\^^ and ako- are the creation and annihilation op- 
erators of a photon with wave vector k and polarization 
a. The electric-field operator is expanded in terms of the 
creation and annihilation operators as 



E(r)=.^ 



2^ 



(2) 
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where Cko- is the polarization vector of photons. The 
internal states of two-level atoms are described by the 
Hamiltonian 



Hf,\ — 



(3) 



where N is the number of atoms, Tiloa is the transition 
energy of the internal levels, and the operator cTj^ is the 
Pauli operator which takes on the eigenvalue 1 (or —1) 
when the jth atom is in the excited (or ground) inter- 
nal level. The Hamiltonian corresponding to the kinetic 
degrees of freedom of the atoms is given by 



ak 



N 

E 



V 

2m ^ 



f:^.(R,-R,), 

(4) 



where m is the mass of the atom and Vt(R) is the trap 
potential. The last term in Eq. (^ describes the Fermi's 
contact type of interaction between hard-core atoms, 
which is characterized by the s-wave scattering length 
a. 

We assume the electric-dipole interaction between the 
atoms and the photon field which is described by the 
Hamiltonian 



N 



Hcd - - ^ • E(Rj). 



j=i 



The dipole operator denoted as D,- has the form 



d(7j+ -I- d*(Tj_, 



(5) 



(6) 



where d is the electric-dipole matrix element, and (Tj+ 
and (7j_ are the raising and lowering operators of the 
internal state of the jth atom. 

The Hamiltonian for the atom-laser interaction is given 

by 

N 



-f£(eL •d)*e*("'^*-'^^•^^)CTJ 



(7) 



where £l, cl, ^^l, and k/, denote the amplitude, the po- 
larization vector, the frequency, and the wave vector of 
the laser field, respectively. Since the optical processes 
of the internal states of the atoms are much faster than 
those of the kinetic motion, we may adiabatically elim- 
inate the dynamics of the internal state. When the de- 
tuning 5 = ujl ^ is much larger than both the atom- 
laser field coupling and the collective linewidth F of the 
atoms 0,p5|], i.e.. 



\E,d\ 



\5\ » r. 



(8) 



the probability that the atoms are in the excited state 
|e ) is much smaller than that in the ground state \g )■ 
In this case, the atomic internal state may be assumed 
to be in the ground state of the Hamiltonian Hg\ H- , 
which is given, up to the first order in jfidj/M, by 



N 

n 



^i(wi,t-ki,-Rj)/2 



9)j + 



\ £lA\ 



-i(wi,t-ki,-Rj)/2| 



(9) 

Taking the expectation value of H^d in Eq- @ with re- 
spect to the atomic internal state (^, we obtain the ef- 
fective Hamiltonian describing the interaction between 
photons and kinetic motion of the atoms via the atomic 
internal levels as 



N 



hU - E 



I^^Ldl 



E(R,). 



(10) 



In the following discussions we will employ the effective 
Hamiltonian HcS = Hf + Hak + H[,^ to describe the dy- 
namics of photons and kinetic motion of the atoms. 

The equation of motion for the electric-field operator 
in the Heisenberg representation is obtained from the 
Hamiltonian H^s as 
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E(r,t) = -[iJ^,(t),V2A(r,t)], (11) 



where the vector potential is given by 



A(r,t)=^ 



k(T 



e' ■''eko-ako-(0 



2eouJkV 
=k.«L(i) ■ 

The equation of motion ([ll| ) can be solved, giving 
E(r,<) =Evac(r,t)+Es(r,i), 



(12) 



(13) 



where Evac(r,t) — en^f*E(r)e r-^'* is the vacuum field 
and 



^'^'''*> - ~ Ann J "^"^ \r-r'\ 
xS{t-t' -\r-r'\/c) 



(14) 



describes the source field scattered from the atoms. Eval- 
uating the commutator and carrying out the integral with 
respect to t' , the positive frequency part of the source 
field (14) becomes 



E(+)(r,<) = - 



1 |gLd| 



dv'- 



|r-r'| 

N 



^ k X (d X k)e-^''-'-' ^ e^(k^+k)-R,(t.)^ 

k j=l 

(15) 
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where tr 
Kjitr) = ei 



t — \r — r'l/c and 
-E^etft.. Since jfidl/M < 1, it is 
sufficient to keep only the terms first order in |£Ld|/?i5, 
which amounts to disregarding the multiple scattering 
and the dipole-dipole interaction. Consequently, we can 
make an approximation 

KjiU) ~ e^^'"<*Rje~^^--'^* = Ilj(t), (16) 

where we ignore the time |r — r'|/c it takes photons to 
propagate from the atoms to the detector as it is much 
smaller than the time scale of the kinetic motion of the 
atoms. The integrand of the source field (|l5|) contributes 
only when r' falls within the spread of the atomic wave 
function. Therefore, if we choose the origin of the coordi- 
nate system at the center of the trap, the distance |r| from 
the atom to the detector is much larger than |r'|, and 
hence — t — (r/c) (l — r • r'/r^) and l/|r — r'| ~ 1/r. 
We thus obtain 



E(+)(r,t) 



ineorV Sh 



-iujL i 



N 



I dr' ^ k X (d X k) 



^^£^n X (d X n)e-^"-*+^'=-'^ 

47reor 5h 

N 

X ^ gi(kt-fctn)-Rj(t) 



N 



F(r,t)^e'^''-^^W, 



(17) 



where n = r/|r| and Ak = k^ — /clii. The expression of 
the source- field operator ( p7| ) consists of the factor F(r, t) 
which depends on the spatial configuration of the system 
and the operator that transfers a recoil momentum ?iAk 
to the atoms. 



III. FIRST AND SECOND-ORDER COHERENCE 
OF PHOTONS SCATTERED FROM A 
NONINTERACTING CONDENSATE 

In this section we consider an ideal Bose gas, where 
the atom-atom interaction is absent. We write the eigen- 
functions of the Hamiltonian i/ak as /„(R), where n is 
an index for the eigenstates with n — Q referring to the 
ground state. 

It is convenient to discuss the problem in the second- 
quantized formalism. We expand the atomic field opera- 
tor as 



i^(R) = ^/„(R)S„ 



(18) 



where 6„ annihilates an atom in the nth eigenstate, and 
satisfies the Bose commutation relation = (5„„'. 

The Hamiltonian ffak can then be rewritten as 



(19) 



where hujn is the nth eigenenergy. The second-quantized 
form of the operator e*'^'''^^ '-*^ in Eq. ( p7| ) is given 

by 

i»(Ak,i) = ydRV>t(R, t)e'^^-'^i}{K, t) 

= ^(n'|e'^'^-^|n)6li„e'('""'-'"")*, (20) 



where V'(r,i) = 
and the bracket {n'\e'^^^''^\n 



^(R)e 



is a shorthand notation of 
/ dR/*,(R)e'^'^ *^/„(R). The source-field operator (|l|) 
is expressed as M+\r,t) = F{r,t)f>{Ak,t). 

At zero temperature all N atoms are Bose-Einstein 
condensed in the ground-state level of the trap poten- 
tial. We assume an isotropic harmonic trap 14 (R) = 
muo^E? /2 for simplicity. The expectation value of the 
electric field at the position r becomes 



(EW(r,t)) =F(r,i)iV(0|e'^''-^|0) 
= F(r,^)iVe-l'^^''l'/^ 



(21) 
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where d = {U/mui) denotes the characteristic length 
scale of the trap. The exponential factor in Eq. (|2l] ) is a 
decreasing function of the scattering angle. 

Suppose that we detect the number of scattered pho- 
tons in a configuration schematically illustrated in Fig. 
The position of the ith detector is denoted as and the 
scattering angle as measured from the direction of laser 
propagation is denoted as 9i. The number of photons 
counted by the detector is proportional to the normally- 
ordered intensity of the field 

/(r,i)^(E(-)(r,t).E(+Hr,0)- (22) 
Substituting Eq. ( [l7| ) into this gives 

/(r, t) = |F(r, t)\^N [l + (A^ - l)e-l'*'^k| V2I _ (23) 

The first term on the right-hand side is proportional to 
the number of atoms N and the angular dependence lies 
solely in |F(r,f)p. The second term is proportional to 
N{N — 1) and contains a forward-scattering factor of 
g-|dAk| /2_ jj-^ ia,ct, in experiments of Ref. ||], a forward 
scattering that satisfies |dAk| < 1 was observed at the 
onset of BEC. The degree of first-order coherence is de- 
fined by pi 



gW(ri,ii;r2,t2) 



(E(-)(ri,ti)-E(+)(r2,t2)) 
V/(ri,ti)/(r2,i2) 



(24) 
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FIG. 1. Schematic of a scattering experiment. The 
trapped Bose-Einstein condensate (BEC) is irradiated by a 
weak and far off-resonant laser. The scattered photons are 
detected by the detectors set around the BEC. The positions 
of the detectors are denoted as and the scattering angles 
as Oi {i = 1,2). 

Substituting Eq. ([l7| ) into this gives 
|5(i)(ri,ii;r2,i2)| 

g-|dAki|V4-|dAk2|V4 



d^Aki-Ak2e-'"<*i-'2)/2 _|_ — 1 



[1 + (iV - l)e-MAki|V2] 2 ^1 ^_ _ i)g-|dAk2|V2] 2 

(25) 

where Aki = — A;Li"i/|i"i| (i = 1 and 2). When 
the number of atoms is sufficiently large to satisfy 
jYg-|dAk.|V2 > 1 for i = 1 and 2, the degree of the first- 
order coherence becomes |(jf^^^(ri, ii; r2, i2)| — 1- The 
scattered photon field is therefore first-order coherent 
when the number of atoms in the condensate is suffi- 
ciently large. On the other hand, when TVe"''''^*''! ^ 1 
for both i = 1 and 2, Eq. (ES) reduces to 



5(^)(ri,ii;r2,t2)| ^ exp I -^|Akie 



-iuj(ti —t2) 



Ako 



(26) 

When ti = t2, the coherence ( p6| ) is lost exponentially 
with increasing the difference in the scattering angles 
|Aki — Ak2|. This is because two photon paths cor- 
responding to atoms kicked in different directions do not 
interfere with each other. The coherence is periodic with 
respect to the time interval ti — t2 with the period of 
27TUJ~^. We note that the same expression as Eq. ( |2^ ) 
can be obtained by substituting iV = 1 in Eq. (p5|), and 
therefore the correlation described above cannot be as- 
cribed to BEC. 

Figure || shows the degree of first-order coherence for 



t2, \g'^^Hri,t;r2,t)\, where N 



1/xm, 



and 27r/A;i — 589nm (the 3s-3p transition of ■'Na) are 



FIG. 2. The angular dependence of the degree of 
first-order coherence g^-^\ri,t;r2,t) of photons scattered 
from 10^ BEC atoms. The wavelength of the laser field 
2-K/kL ~ 589nm (the 3s-3p transition of ^''Na) and the char- 
acteristic length of the trap d — Ifim are assumed. The wave 
vectors ki and k2 are assumed to be in the same plane, and 
9i and 62 denote their angles as measured from k^. 

assumed. The wave vector ki and k2 are assumed to 
be in the same plane, and 61 and 62 denote the angles 
of ki and k2 as measured from ki (see Fig. ||). One 
can see that g*^^-' ~ 1 for jdAkj < 1, and therefore, the 
photons scattered in the forward directions have first- 
order coherence. There are peaks along 9i — O2, which 
can be seen from Eq. (|26|). 

The degree of second-order coherence p4] is defined by 



g(2)(ri,ii;r2,i2) 



{eL \ri,h)E}j 'ir2,t2)E}^>ir2,t2)E'J'iri,h)) 



^(-) 



a..f3—x.y,z 



(27) 



/(ri,ii)/(r2,i2) 

Substituting Eq. ( p7| ) into this yields 
g(')(ri,ii;r2,i2) = 

[n{N - 1){N - 2){N - 3)e-MAki|V2-|dAk2|V2 



+N{N - 1){N ~2) 
xRe(e-'^'^''i- 



2e 



(iAkip/2-|ciAk2p/2 



Ak2e''^<*2-'i)/2 _|_ g!i^Aki-Ak2e" 



*'(*2 -'1) ^2^ 



_g-|rfAki|V2 _^g-|rfAk2|V2 



+N{N - 1) 



1 + e" 



|Aki/2-Ak2e*"('2-'i)|2 



_g-rf2|Aki+Ak2e""'*2-tl)|2/2 ^ 2e"l'*^'^il^/2 

_2g-|rfAk2|V2 cos(d2Aki • Ak2 sinw(t2 - h)) + ivj 



/ {tV^ [1 + (iV - i)e-l'''^ki|V2j Ji + (iv - l)e 



dAk2p/2 



(28) 
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If the number of atoms is sufRciently large to satisfy 
jYg-|dAk,|V2 ^ I foj. botij j = 1 and 2, the degree of 
second-order coherence becomes g^^\ri,ti;r2,t2) ~ 1 as 
well as \g^^\ri,ti;r2,t2)\ — 1- In general, rth-order co- 
herence becomes unity for all r, and therefore, in 
this large TV limit the scattered photons are coherent in 
all orders. On the other hand, if ATe-M^^ilVz ^ i Jq^. 
both i = 1 and 2, Eq. ( p8| ) reduces to 

^_g-d"|Aki+Ak2e"*'<'2-*i)|2/2^ ^29) 

When ti — t2, Eq. ( p9| ) becomes g^'^^ ~ 2 for either 
Aki ~ Ak2 or Aki ~ — Ak2. This means that a pair of 
photons tends to be scattered either in the same direc- 
tion or in the symmetrical directions with respect to laser 
propagation. These correlation in scattered photons can 
be interpreted as follows. In the case of Aki — Ak2, the 
first photon kicks an atom from the condensate to the 
state e*^''i'^/o(R), and so does the subsequent photon, 
where /o(R) is the ground state of the trap potential. 
Consequently, the Bose-enhancement factor of this pro- 
cess is 

((7V-2)o,2AkJ6ik.^ofoLi^o|iVo,OAkJ ^ ^2N{N-l), 

(30) 

where bl = E„("|e*'^|0>&]; and |mk) = b]^/V^.\0). 
Therefore the probability that this process occurs is en- 
hanced by a factor of two. In the case of Aki — — Ak2, 
there is a process in which the subsequent photon kicks 
back the atom excited by the first photon to the con- 
densate, in addition to the process in which two atoms 
are excited to the states e*^'^i -^/o(R) and e*'^''^-^/o(R). 
The Bose-enhancement factor of each process is 

(iVo,0AkJfeS6Aki6Lki^o|A^o,0Aki) - N (31) 

and 

((iV-2)o,lAki,lAk.|&SfoAki ^Aki ^0 1 ^0 , Aki , Ak2 ) 

~ ^NiN-1), (32) 

and hence the probability is enhanced by a factor of two 
for large N. Thus the above correlations are due to the 
quantum-statistical effect of bosons. In fact, the degree of 
second-order coherence reduces to g'^^ = 1 when iV = 1, 
and a photon pair scattered from a single atom shows no 
correlation. We note that the presence of the condensate 
is required for the above discussion. The process in which 
the condensate participates is accompanied by the Bose- 
enhancement factor of ^/N, and the probability that the 
process occurs is N times greater than that of the process 
between other levels. Therefore, scattering of photons by 
the condensate can clearly be distinguished from that by 
the non-condensed atoms. When the condensate is not 




FIG. 3. The angular dependence of the degree of 
second-order coherence g^^\r-i,t;r2,t) of photons scattered 
from 10'' BEG atoms. The wavelength of the laser field 
2Tv/kL ~ 589nm (the 3s-3p transition of ^'^Na) and the char- 
acteristic length of the trap d — Ifim are assumed. The wave 
vectors ki and k2 are assumed to be in the same plane, and 
di and 02 denote their angles as measured from k^. 

present and the atoms are thermally distributed, occu- 
pation numbers of any states are smaller than one, and 
Bose-enhancement effects cannot be observed. 

The degree of second-order coherence (^8|) for ii — t2 
is shown in Fig. ^, where N = 10^ and d — 1/im, and ki 
and k2 are assumed to be on the same plane. We find 
that g^^^ ~ 1 when the scattering angles are small, and 
gi^) grows with increasing 9i. While the peak of g'-^-' at 
the same scattering angle 6i — 62 is g'^^^ ~ 2 for large 
9i, the peak at the opposite angle 9i = —62 first grows 
and then decreases to g'^'^^ ~ 1 as 0i increases. This is 
due to the fact that Aki = — Ak2 is not compatible with 
the energy conservation |ki| = |k2| = (^l/c (the recoil 
energy is much smaller than the photon energy) , and the 
third term in Eq. ( p9|) cannot become unity even for 9i = 
—92- In other words, the energy conservation prohibits 
the second photon from exactly kicking the atom back to 
the condensate. 

Figure ^ shows the time dependence of the degree of 
second-order coherence (p8|). The parameters are the 
same as those in Fig. |3[ and the scattering angle of the 
first photon is fixed at 9i = 7r/12. We find that the 
second-order coherence is periodic in t2 — ti with the pe- 
riod given by 'r:uj~^. This periodic behavior reflects the 
fact that a recoil atom oscillates in the trap. The co- 
herence disappears when the atom is kicked out of the 
condensate, and revives when the atom returns to the 
condensate. We note that g^^^ goes to below unity be- 
fore and after the appearance of the peaks. This indi- 
cates that photons are suppressed to be scattered into 
6*2 — ±01 in these time intervals. This comes from the 
terms proportional to N{N — \){N — 2) in the numera- 
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FIG. 4. The time and angular dependence of the degree 
of second-order coherence g^^' (ri, ti; r2, t2) of photons scat- 
tered from 10'^ BEG atoms. The scattering angle of the first 
photon is fixed at 9\ = it/12. The wavelength of the laser 
field 2-K jkh = 589nm (the 3s-3p transition of ^^Na) and the 
characteristic length of the trap d — 1/im are assumed. The 
wave vectors ki and k2 are assumed to be in the same plane, 
and 6i and 62 denote their angles as measured from k_L. 



tor of Eq. (28), which shows destructive interference of 
probability amplitudes of two-photon processes in which 
one photon kicks an atom and another photon does not 
disturb the atoms. 

The degree of first-order coherence 3*^^^ can be mea- 
sured by superposing the photon field in the directions 
of ki and k2 with a delay time t2 — ti in the path of the 
ki photon. Taking a delay of order lo^^ ^ 1 ms while 
retaining the coherence might be feasible by reducing the 
speed of light using an ultracold atomic gas . On the 
other hand, a measurement of g'-^-' does not require the 
coherent delay, since it measures a correlation between 
detected photons. 



IV. EFFECTS OF ATOM-ATOM INTERACTIONS 

In this section, we take into account the atom-atom 
interaction on the coherence properties of photons scat- 
tered from BEG within the Bogoliubov approxima- 
tion ||. 

The second-quantized form of the kinetic part of the 
Hamiltonian is given by 



dR?A^(R)V't(R)7/.(R)v;(R; 



V(R) 



Uo 
2 



(33) 



where V'(R-) is the field operator of the atoms ( |18| ) and 
Ua = 4:TTh^a/m, a being the s-wave scattering length. In 



the Bogoliubov approximation, the atomic field operator 
is divided into a mean field part and fluctuations from it 
as 

(34) 



where V'g(R) is the normalized single-particle wave func- 
tion which obeys the Gross-Pitaevskii (GP) equation |23] 



2m 



V^ + Vt{R) + UoN\i;{R)\' 



i;{R) = /.^(R) 



(35) 



where /i is determined so that the solution of Eq. (|35| ) 
satisfies the normalization condition J (iR|?/'(R)p — 1. 
Equation (|3^ ) guarantees that if we substitute Eq. (jsj) 
into Eq. (|33|), terms linear in S^p or 6^p^ in Eq. (|3^) van- 
ish identically. Because the interaction between atoms 
is weak, we may expect only small deviations from the 
mean field at zero temperature, so that we may keep only 
terms up to quadratic in Sip and 6ip^: 



K c^N J dRiP*g{R) 



2m 



V's(R) 



Vt(R) 



' 2m 



VAR) 



dR 

+2UoN\il;{R)\^jS'tP{R) 
+UoNi^f{R)diP^R) + UoN^l{R)SiP^^{R) 



(36) 



To diagonalize the quadratic terms, we write Sip{R) as 



5V'(R) = E["^(i^)'^A + «I(R)/3l 



(37) 



Requiring the coefficients of the terms and p]^ to van- 
ish, we find that ua(R) and wa(R) should satisfy 



— + VtiR) -fi + 2UoN\iPg{R)\'' 
Zm 

+UoN^g{RfvxiR) = hujxuxiR), 



uxiR) 



(38a) 



2m 



V^ + Vt{R)-^l + 2UoN\iJg{R)\' 



vx{R) 



+UoNr'{R)ux{R) = -hLOxvxiR). 



(38b) 



The quadratic part of the Hamiltonian ( p6D can then be 
diagonalized as 



(39) 



where f3x are the creation and annihilation operators 
of the Bogoliubov quasiparticles and huj\ denotes the en- 
ergy of elementary excitations. The constant terms are 
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irrelevant to later discussions and are therefore omitted 
in Eq. 

The second quantized form of the op- 
erator 6*'^''''^^^*'' that gives the momentum to an 
atom reads in the Bogoliubov approximation as 



I?B(Ak,t)= / dRi/)T(R,t)e'^'^-^7/;(R,t), 



(40) 



where V'(R,i) = e^'^^iJ;(R)e~i^\ Because |ma(R-)| > 
|ua(R)| pel we may neglect vx{R) in Eq. (0. The 
operator (Hg) then becomes 



Ps(Ak, t) 



„jAkR| 



u\)(3\e 



-iuJxt 



4^ 




+ {u\\e 



lAkR 



\ipg) Px^^'^^^l FIG. 5. The angular dependence of the degree of 

first-order coherence g'^' (ri , t ; r2 , t) of photons scattered from 
|giAk-R|y^^^t^^^gi(a;A-c:jf,)t^ (41) 10^ interacting BEG atoms. The s-wave scattering length is 



where we use the shorthand notation as in Eq. (20) 



In the Bogoliubov approximation, the orthonormality 
condition is expressed as {u\\u^) — {v\\v^) = (5a/j, and 
the completeness relation is given by J2xiW>^){^^\ ~ 
\v\){v\\) = 1 — \^g){^pg\. In the present approxima- 
tion of neglecting v\, they reduce to {ux\u^) = and 



Ea \u\){ux\ + \i'g){^gl 

Using the operator 



= 1, respectively. 
tl|), the source-field operator is 

expressed as Er^(r,i) = F(r, t)X'B(Ak, t). The expec- 
tation values of the electric field and the intensity are 
calculated to be 



(E(+)(r,i)) =F(r,t)iV(V3 



lAkR 



1^: 



3/' 



I{r,t) = \F{r,t)\'N{l + N\{^g\ 



iAk-R 



1^. 



(42) 



(43) 



These expressions are different from the noninteracting 
counterparts (|l]) and ( p3| ) in that the expectation value 
with respect to the ground state of the noninteracting 
|0) is replaced by that of the ground state 



atoms 



iAk-R 



of the GP equation {iPg\e'^^-'^\iJjg) , and - 1 by N. 
The latter arises from the assumption of the Bogoliubov 
theory that the ground state is a coherent state rather 
than a Fock state. Although this is an artifact, the error 
is of order 1/N and negligible for A'' ^ 1. The former 
reflects the fact that the ground-state wave function of 
interacting atoms expands for the repulsive interaction 
a > 0, and contracts for the attractive interaction a < 0, 
compared with that of noninteracting atoms. 

The degree of first-order coherence is calculated to be 



g(i)(ri,ii;r2,i2)| - 

^g-»Aki.R^^(^^ _^2)g»Ak2-R^ 



+ N{e 



-iAki-R\ 



iAk2-R\ 



assumed to be a = 2.75 nm. The wavelength of the laser field 
2-K/kL ~ 589nm (the 3s-3p transition of ^''Na) and the char- 
acteristic length of the trap d = l/^m are assumed. The wave 
vectors ki and k2 are assumed to be in the same plane, and 
6*1 and 62 denote their angles as measured from k^. 



where Usit) = Ea \ux)e~''^'\uxl ^..^ \- - ig 
expectation value with respect to the ground state of the 
GP equation (^5|) . To evaluate the expectation value that 
includes J7_b, we need to find ma which we obtain by nu- 
merically diagonalizing Eqs. (^8^) and (|3^) using the 
method discussed in Ref. |2 7 [ . Figure || shows the degree 
of first-order coherence (|44|) for ti — t2. The parameters 
are taken to be iV = 10 ,d ~ 1/im, and a = 2.75 nm. 
The ground state of the GP equation (p5|) is obtained nu- 
merically |2^ . The sharp dips in Fig. |5| correspond to the 
points at which the sign of g^^^ changes (note that Fig. ^ 
displays the absolute value of 5*-^^), while the curves in 
the noninteracting case are smooth (sec Fig. ||). To un- 
derstand the oscillatory behaviors of g''^^\ we compare 
in Fig. ^ the profile of the interacting ground-state wave 
function with iV = lO'^ with the corresponding nonin- 
teracting one. As we can see from the figure the den- 
sity profile near the origin becomes flat due to the re- 
pulsive interaction which gives rise to the oscillatory be- 
havior of e*^"^'^. The validity of this interpretation can 
be confirmed analytically in the limit of strong interac- 
tion Na/d 3> 1. In this limit the kinetic energy may be 
neglected in comparison with the trap potential and the 
interaction energy (Thomas-Fermi approximation), and 
the wave function is given by 



and ' 



denotes the 



V'TF(i?) = 



/ 15 



(45) 



(l + ^Ke^'^ki.R^^p)^ (l + iV|(e 



,iAk2-R\ 



,P)' where Rq = d{15Na/d)^^^ In this Thomas-Fermi 

(44) limit, the expectation value becomes 
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0.3 
0.2 
0.1 




\< noninteracting 



interacting 








r/d 



FIG. 6. The profiles of tlie ground-state wave functions. 
Tlie solid curve shows the wave function of interacting 10"^ 
atoms which is obtained by solving the GP equation (^5[), 
and the dashed curve shows the noninteracting one. 



15 



[(6-2i?,o|Ak|)sin i?o|Ak| 
(46) 



2|i?oAk|5 

-6i?o|Ak|cos i?o|Ak|], 



which shows oscillations of g'-^K In contrast the Gaussian 
wave function gives a monotonically decreasing nonzero 
value (e'^"^-^) = e"!''^'^!'/*. 

The degree of second-order coherence is calculated to 

be 



4i/„iAki-R\ / lAka-Rx |2 
/ffW /g\ 

-iAki-R\ /„-iAk2-R\ 



iV-^ 2Re(e" 



X e 



iAk2R 



-h2Re(e-'^'^i-^)g(e'^'^^-^) 



„-iAk2-R 



UB{t2 ~ ti)e'^''^-^) 



,iAki-R\ |2 



,jAk2R\ 



/Si /g\ 



/lAki-R^ |2 



+N^[l + \{e''^''--'^UB{t2-h)e' 

^|^g-,Ak2.Rjj^(^^„^^)g»Ak,.R^^|2 

_^2|(e»Aki.R^^|2 ^ 2Re{e-'^^'-^)g 



+N 



(1 + N\{e"^^'-'^)g\^) (1 + N\{e"^^'-'^)g\^) 



(47) 



Figure ^ shows the degree of second-order coherence ( p7| ) 
for — t2, where the parameters are the same as in 
Fig. ^. One can see the complicated structure that arises 




FIG. 7. The angular dependence of the degree of sec- 
ond-order coherence g^^\ri,t;r2,t) of photons scattered from 
10^ interacting BEG atoms. The s-wave scattering length is 
assumed to be a = 2.75 nm. The wavelength of a laser field 
27r/fc_L = 589nm (the 3s-3p transition of ^''Na) and the char- 
acteristic length of the trap d = Ifim are assumed. The wave 
vectors ki and k2 are assumed to be in the same plane, and 
6*1 and 62 denote their angles as measured from k^. 



from the oscillations of (e*'^''^ '^)^. Figure || shows the de- 
pendence on the time interval t2 — ti of (7^^^(rijti; r2, t2), 
where the parameters are the same as in Fig. H, and the 
scattering angle of the first photon 9i is fixed at 7r/12. We 
find that the peaks appear periodically with their heights 
monotonically decreasing. The decay of the peaks can be 
understood as follows. The atom is kicked out of the con- 
densate upon receipt of recoil momentum, and oscillates 
in the trap potential while encountering the condensate 
periodically. If the atom-atom interaction is not present, 
the Gaussian wave function of the scattered atom retains 
its shape during the oscillations, and the peaks revive 
completely as shown in Fig. ^ On the other hand, when 
the atom- atom interaction is present, the wave function 
undergoes distortion by colliding the condensate. As a 
result, the atom cannot be kicked back completely to the 
condensate by scattering a photon, and therefore the cor- 
relation deteriorates. 



V. CONCLUSIONS 

We have studied the coherence properties of photons 
scattered from the BEG irradiated by a weak and far- 
off-resonant laser field, and have shown that the degree 
of second-order coherence g'--^"^ exhibits unique features 
of trapped BEG. That is, the probability of two photons 
being scattered either in the same direction or in the sym- 
metrical direction with respect to laser propagation is en- 
hanced, and these correlations appear periodically with 
the period of atomic motion in the trapping potential. 
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FIG. 8. The time and angular dependence of the degree of 
second-order coherence g''^' (ri, ii; r2, i2) of photons scattered 
from 10^ interacting BEC atoms. The scattering angle of the 
first photon is fixed at 6i = it/12. The s-wave scattering 
length is assumed to be a = 2.75 nm. The wavelength of the 
laser field 2-R/kL = 589nm (the 3s-3p transition of ^''Na) and 
the characteristic length of the trap d — l/im are assumed. 
The wave vectors ki and k2 are assumed to be in the same 
plane, and 6i and 02 denote their angles as measured from 

These correlations in space and time are manifestations 
of the combined effect of bosonic stimulation and quan- 
tized kinetic motion in the trap. We have also taken into 
account effects of the atom-atom interaction, and have 
confirmed that such correlations can be observed even 
in the presence of the atom-atom interaction with ap- 
pearance of an additional structure (sharp dips in Fig. ||) 
resulting from the repulsive nature of the interaction. 
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